Genereer een gemeentekaart van Nederland

De data is online beschikbaar. O.a. (hier)[https://data.overheid.nl/data/dataset/grenzen-van-alle-nederlandse-gemeenten-en-provincies/resource/d5570855-ceb7-458b-ba93-61e5ec332649]. Plaats alle shapefiles in dezelfde folder als je R script en lees de data in.

library(rgdal)
gemeentenShp <- readOGR(".","gemeenten")
## OGR data source with driver: ESRI Shapefile 
## Source: ".", layer: "gemeenten"
## with 393 features
## It has 6 fields

Converteer de coordinaten naar een voor R bruikbaar format.

gemeentenShp <- spTransform(gemeentenShp, CRS("+proj=longlat +datum=WGS84"))

Zet de lijst om naar een data frame.

library(broom)
gemeentenDF <- tidy(gemeentenShp)

Bouw de kaart.

library(ggplot2)
library(ggmap)
nl<-get_map("Netherlands",zoom=8)
ggmap(nl) + 
  geom_polygon(aes(x=long, y=lat, group=group),
               fill='grey',
               size=.2,color='green',
               data=gemeentenDF,
               alpha=0)

Source:

D. Kahle and H. Wickham. ggmap: Spatial Visualization with ggplot2. The R Journal, 5(1), 144-161. URL http://journal.r-project.org/archive/2013-1/kahle-wickham.pdf